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Abstract 

According to a proposal by Arnold and Yaffe, the non-perturbative g 2 T-contribution to 
the Debye mass in the deconfmed QCD plasma phase can be determined from a single 
Wilson line operator in the three-dimensional pure SU(3) gauge theory. We extend a 
previous SU(2) measurement of this quantity to the physical SU(3) case. We find a 
numerical coefficient which is more accurate and smaller than that obtained previously 
with another method, but still very large compared with the naive expectation: the 
correction is larger than the leading term up to T ~ 10 7 T C , corresponding to g 2 ~ 0.4. 
At moderate temperatures T ~ 2T C , a consistent picture emerges where the Debye 
mass is ~ 6T, the lightest gauge invariant screening mass in the system is ~ 3T, 
and the purely magnetic operators couple dominantly to a scale ~ 6T. Electric (~ gT) 
and magnetic (~ g 2 T) scales are therefore strongly overlapping close to the phase 
transition, and the colour-electric fields play an essential role in the dynamics. 
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1 Introduction 



An important concept in the description of the deconfined finite temperature QCD 
plasma is the Debye mass mo- It is the inverse of the spatial screening length felt 
by colour-electric excitations, in analogy with Debye screening in abelian plasmas, 
which is felt by electric fields but not by magnetic fields. The leading order result for 
ran is perturbatively computable 0, but the next-to- leading order result is already 
non-perturbative 0, making the determination of mp (and, in fact, already its precise 
definition) a challenging theoretical problem. The numerical value of has a number 
of applications, in particular in the phenomenology for quark-gluon plasma formation 
(for a review, see 0). It is also interesting to note that, even though mo is originally a 
strictly static quantity, it might have significance even for some dynamical observables, 
such as the rates of anomalous processes in the plasma 

A gauge invariant, non-perturbative definition for wid has been given by Arnold and 
Yaffe ||. This definition is based on a (Euclidian) time reflection operation, denoted 
also by TZ Z , which reverses the sign of the temporal component of the gauge field: 
1Z z Aq — > —A% (for a more detailed discussion, see ||). The Debye mass is defined as the 
smallest among the exponential fall-offs that are found for gauge invariant operators 
odd under this symmetry. A further observation is that due to finite temperature 
dimensional reduction 0-0, the measurement of mj) in numerical simulations can 
be carried out in a three-dimensional (3d) SU(N) + adjoint Higgs theory. Lattice 
measurements in this theory were performed, both for SU(2) and SU(3), in ||10||P|. 

However, there is a problem in making precise simulations in the 3d SU(3) + adjoint 
Higgs theory. The reason is that in the regime of asymptotically small couplings, 
there is a scale hierarchy between the leading order result tub ~ gT and the non- 
perturbative correction Sntp ~ g 2 T, which one is trying to determine (the "magnetic" 
degrees of freedom are also associated with the scale g 2 T). As always, a scale hierarchy 
is difficult to accommodate on a single lattice of presently available sizes, and thus 
makes the approach to the continuum limit quite problematic in practical simulations. 

Fortunately, there is a way out: parametrically (i.e., as a power series in g), even the 
field Aq can be integrated out perturbatively, leaving a pure 3d SU(3) gauge theory. 
The operator odd in Aq goes over to a non-local but gauge invariant Wilson line in the 
pure SU(3) gauge theory 0. This allows to determine the coefficient of g 2 T. 

For the SU(2) case, this Wilson line operator was studied with lattice simulations 
in ]TJ| . We found a result for the g 2 T term which is more precise and smaller than that 



found in flOR . The purpose of the present paper is to extend the SU(2) measurement 
to the physical SU(3) case. 

3 Progrcss towards a lattice measurement of mo directly with 4d lattice simulations has been 
reported in jll]] (see below). Another gauge invariant definition within the 3d SU(3)+adjoint Higgs 
model was suggested in Q. Other definitions of ma, based on gauge fixing, have recently been 
pursued in |l3fl. 
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2 The operator in the continuum 



Let us start by briefly reviewing the basic idea of the approach in [H. Consider gauge 
invariant operators 0(x) odd under A — > —A (i.e., in the notation of f5|, with 1Z Z = 
— 1). In practice, we consider 0(x) containing only one power of Aq, since otherwise 
(for Aq , n > 3) the leading behaviour for small g is expected to be of the form n ■ mo, 
where mo = {N/3 + A^/6) 1//2 gT is the leading order Debye mass (Nf is the number 
of flavours). Imagine now doing the path integral over A before doing that over Aj. 
Then, 

JVA exp(-S Ao ) 

where O' = dO/dA , 



S 



An 



J d 3 x [-(AA) a (AA) a + -m^AgAg] , (2) 



and G ' (x,y) is the tree-level Green's function in the background of Af. 



~6/„ .a exp(- m o|z-2/|) T 



G»*(s,y) = " v U| ^ a& (x,y), (3) 
47r|x — y| 

with W ab {x, y) the adjoint Wilson line, 

W ab {x,y) = 2TiT a W lund (x,y)T b [W^ nd {x,y)}^ (4) 

where W iund (x,y) = V exp(ig \ v x dxiA?T a ) . 

Thus, the measurement of the Debye mass can be reduced to the measurement of 

G F (x,y) = (0' a {x)W ab {x,y)0' b {y)), (5) 

in the pure SU(N) gauge theory. The measurement of this operator produces an expo- 
nential decay with a coefficient m, to which the perturbative part m from Eq. (|3|) has 
to be added in order to obtain the physical result. 

Including also loop corrections (i.e., allowing A { to fluctuate at all length scales), 
maintains the form of the answer but replaces the perturbative mo in Eq. @ by an 
effective parameter mf M: 

,eff__ , N 9if,m \ ( 2 m g\ 



m - mf = m + (h^ + d A ) + 0[g^, ^, ...j, (6) 

where g\ = g 2 T, A is some ultraviolet cutoff, and d\ depends on A. In particular, in 
the case of a lattice regularization, A = a -1 and d\ = —(3 In 2 + l)/2 (see [|, |14|). It 
follows that if the physical Debye mass is written as 



m D = m + ^\n^ + c N g 2 3 + <D(g 3 T), (7) 
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then 



m N ( N 
c N = — + — \\n 



where m is the mass measured from the operator in Eq. (j5|), and (3q = 2N/(ag^). 

We thus need to find operators O' such that we get the lowest possible coefficient 
for the exponential falloff in Eq. (|5|). From | 10| , we expect the operators for Eq. ([[]) to 
be O ~ e 3 jkTr A Fjk and suitable variants thereof (see below), leading to 0' a = F^ k in 
Eq. (||). In principle, one could get rid of the specification of the operators altogether by 
measuring the coefficient of the perimeter law of the adjoint Wilson loop ||. However, 
in practice the perimeter law can only be observed on the lattice by including this very 
kind of operators in the basis, due to a very small overlap of the lattice Wilson loop 
with the asymptotic lightest state |T5 



3 A lattice formulation 

On the lattice, we replace Eq. (|5|) with operators of the type 

GfM^V) = (^(F l3 (x)T a )T a \x,y)i:T(F kl (y)T b ) 
where Fy is the naive lattice analogue of the field strength tensor, 

Fij{x) — — — Pijix) — P^jix) , 

Pij is the plaquette, and T ab is the lattice Wilson line: 

T ab (x,y) = 2Tr(T a S{x,y)T h S\x,y)), 

JV-l 

S(x,y) = Y[U 3 (x + n3), y = x + N ■ 3. 



n=0 



(9) 
(10) 

(11) 
(12) 



Here Z7, is a link variable and 3 is the unit vector in the X3-direction. Eq. (^) may be 
rewritten as 



G 



F,ijkl 



2(Tr 



1 



FMSFMS^ --TiF l0 {x)TiF kl {y) 



N 



(13) 



where the latter term vanishes fast close to the continuum limit and is not numerically 
important (for SU(2), it vanishes identically). In our simulations, we use Gf in this 
last form, with the components ij = kl = 12. 

In our practical simulations, each plaquette at the end of the Wilson line was replaced 
by the sum over all four spatial plaquettes of the same orientation sharing the end point 
of the Wilson line, a configuration called the "clover" . This operator is expected to have 
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Figure 1: The local mass values, ma = — ln[Gp(n + 1)/Gf(n)], from typical diago- 
nalized operators at the different /3(j's. The masses have been converted to continuum 
units by multiplying the lattice values ma with /3g/6. 



a better projection onto the ground state (i.e. to be less contaminated by higher spin 
excitations as well as to improve on statistics. We have carried out measurements 
using single plaquettes, as well, and checked that the two operator types give consistent 
results, with the clover indeed providing a much better signal. In order to further 
increase the projection onto the lightest mass, we used smearing and diagonalization, 
The idea is to measure a whole cross correlation matrix within a channel of 



as m 
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given quantum numbers, and then to search for the optimal operator by a variational 
calculation. The basis is obtained from gauge invariant "smeared" operators with 
different numbers of smearing iterations. For details, we refer to 
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We gathered statistics until the statistical errors of the local masses at their plateaux 
were well below 10%, which required between 500 and 4000 measurements, depending 
on the lattice size. The choice of lattice sizes was governed by the experience from 
SU(3) glueball measurements |16|] as well as our previous SU(2) measurements [II 



We have also explicitly checked that finite volume effects are within the statistical 
errors at single values of (3q (see Tables [I], |2|), and then increased the volume with j3c, 
keeping it fixed in physical units. 
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0.397(3) 


1.59(2) 


1.12(2) 


20 


52 3 


0.327(5) 


1.64(3) 


1.13(3) 


25 


64 3 


0.268(2) 


1.68(2) 


1.14(2) 


oo 








1.14(4) 



Table 1: The results for SU(2). The continuum extrapolation is from a linear fit to 
p G > 12, with x 2 /d.o.f. = 0.03 (see Fig. 0). 

4 Numerical results 

Typical local mass values from the correlation functions in the case of SU(3) are shown 
in Fig. |I|, exhibiting rather pronounced plateaux. The mass values quoted in this paper 
have been obtained by fitting an exponential to the optimal diagonalized correlation 
function, and choosing among various fitting ranges those giving the best values for 
X 2 /dof. The lengths of the fitting ranges are between 5 and 11 timeslices, and the 
X 2 -values are between 0.5 - 1.5. We now discuss our results both for SU(2) and SU(3). 

SU(2). Previously we carried out measurements for SU(2) with 0g = 9, 12, 16 fl4|| . 
We have now added values closer to the continuum limit, 0g = 20, 25. The results are 
shown in Table [I] and in Fig. As a continuum extrapolation we obtain c 2 = 1.14(4). 
Note that we see the onset of a quadratic 0(a 2 ) term in c 2 on coarser lattices (the 
same is true also for SU(3)). 

SU(3). The results for SU(3) are shown in Table and in Fig. g As for SU(2), the 
results show could scaling with (3q after the removal of the logarithmic divergence. As 
a final continuum value we cite c 3 = 1.65(6). 



5 Discussion 



Comparison with other measurements in 3d. Let us compare our results with 
those of previous studies. In [|T0|| , c 2 , c 3 were measured using operators O ~ £3jfcTr A Fjf. 
in the 3d SU(N) + adjoint Higgs theory, with the result c 2 = 1.58(20), c 3 = 2.46(15). 
In [l4|, we found the smaller result c 2 = 1.06(4), using the field strength correlator 
in the pure 3d SU(2) theory as in this work. Adding two more lattice spacings to 
SU(2) and evaluating also SU(3), our present results are c 2 =1.14(4), c 3 =1.65(6). 
In conclusion, the values are considerably smaller than in []10| . However, there the 



system had multiple length scales, so it need not be a surprise that the result has some 
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Figure 2: The results for SU(2) (left) and SU(3) (right). The continuum fits from 
Tables p], have also been shown. The ranges of 1/Pg on the x-axes have been chosen 
such that in lattice units, the masses have comparable values for SU(2) and SU(3). 
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Table 2: The results for SU(3). The continuum extrapolation is from a quadratic fit 
to all the Pg~ values, with y 2 /d.o.i. = 1.25. We use a quadratic fit for SU(3) since we 
are slightly farther from the continuum limit than for SU(2), see Fig. |2|. A linear fit to 
Pa > 28 gives 1.65(5), with x7d-o.f. = 0.04. 



systematic uncertainties. Our present results scale quite well with N, a behaviour which 



is also observed for the glueball spectrum in SU(N) theories f!B |. Furthermore, the 
bare numbers vaj g\ extracted from the Wilson lines exhibit the logarithmic divergence 
expected from lattice perturbation theory ||, adding further confidence to them. 

The mass scales in units of temperature. We now proceed to estimate the 
numerical value of the Debye mass in units of T. For definiteness we take Nf = 0, 
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in which case T c /Aj^g pa 1.03(19) for SU(3) |p7| . (Let us stress that the values of c 2 , c 3 
as such do not depend on the fermionic content of the original 4d theory, and thus C3 
applies also to the physical QCD with Nf 7^ 0.) At moderate temperatures of, say, 2T C , 
we then estimate that g\ ~ 2.7T ||, and thus mo ~ 6T. This estimate is of course 
somewhat rough, since terms of order g 3 T etc. have been neglected^. However, the 
results in [TO] suggest that these corrections should come with very small coefficients. 

An important question for the construction of effective theories concerns the degree 
to which the Debye mass is perturbative, and the related issue of how well the electric 
(~ gT) and magnetic (~ g 2 T) sectors are separated, once non-perturbative coefficients 
are taken into account. Even at very large temperatures, say, T ~ 10 4 T C with a weaker 
coupling g\ pa 0.64T, the correction (~ LIT) is still larger than the leading term 
(~ 0.8T). Both are of the same order (~ 0.7T) only at T ~ 10 7 T C . 

It is also interesting to compare with the masses in the magnetic sector, especially 
at low temperatures T ~ 2T C . There the ordering of states according to the naive 
expectations is reversed. A purely magnetic excitation, the lightest glueball, has a 
- 6T 



mass 



16 



making it as heavy as mo- However, there is a lighter 
excitation with the same quantum numbers which then determines, e.g., the decay of 
the Polyakov loop. In the 3d SU(3) + adjoint Higgs theory, it is represented by the 
mass of the bound state AqAq, whose numerical value is about ~ 3T f|, 

Comparison with measurements in 4d. All these numbers are in good qual- 
itative agreement with those found with direct 4d simulations in 

m(Af + ] 



11 



of 11 



In Table 5 
channel fin 



2.60(4)T corresponds to the lightest mass in the 
the language of the 3d theory, to ui^AqAq) ~ 3T], while the 1Z Z = — 1 || |TT]] mass 
m(Ai + ) pa 6.3(2)T should correspond to the lightest mass in the electric sector [in the 
language of the 3d theory, to as defined here, ~ 6T]. Since the 4d measurement of 
mo contains all orders in g this agreement would indeed suggest corrections of order 
g 3 T and higher to be small. 

the lightest masses 
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Incidentally, we believe that, contrary to the suggestion in 
measured in [[Ll[] do not correspond to the glueball masses measured in but to scalar 
states involving Aq: the screening masses corresponding to static glueball correlators 
are heavier. Thus, close to T c , the naive expectation that Aq is heavy and decouples 
is completely unjustified (see also || [| |2"U||). This should not be a surprise, since 
the actual transition is assumed to be driven by the Z(N) symmetry related to the 
Polyakov line and Aq. At the same time, it is interesting to note that the higher lying 



glueball spectrum itself is expected to be relatively little affected by Aq [19 



Conclusions. In conclusion, we have determined the non-perturbative contribution 
of the soft modes ~ g 2 T to the Debye mass in QCD and found it to be larger than the 

4 Note also that in principle, some other 1Z Z = — 1 operator such as Tr Aq could have a lighter mass 
in this region. 
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"leading" perturbative contributions for all reasonable temperatures. The approximate 
agreement between the values for the Debye mass as determined in this simulation and 
with the 4d measurement, further suggests that the 0(g 3 T) and higher corrections 
in the Debye mass are small even at relatively small temperatures, as already argued 
jT0|. The consistent picture emerging from our results, as well as from previous 



m 



analytical work and simulations [^, [| [TD|, OT, is that a dimensionally reduced theory 
gives a reliable description of correlation functions down to quite moderate tempera- 
tures T ~ 2T C , when A Q is kept in the action. However, at such temperatures A may 
not be integrated out, but is in fact an essential constituent in the lightest physical 
degrees of freedom. 
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